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TECHNICAL MEMORANDUM X- 64838 


A STOKES' VELOCITY PHOTOGRAPHIC METHOD FOR MEASURING 
THE SIZE DISTRIBUTION OF AEROSOLS 

INTRODUCTION AND THEORY 


The purpose of this research was to provide a (absolute) method for 
measuring aerosol size distributions that was independent of the light scattering 
properties, which form the basis of most commercial optical particle sizers, 
of the particles. Such a method could therefore be useful in calibrating and 
testing the performance of an optical scattering sizer that has been designed 
and built at Wayne State University 1 in an effort to provide standard sources 
of known size particles for calibrating LDV systems. 2 

The basis of the method is (die Stokes’ law [ ll , which relates the ter- 
minal velocity (v ) of a sphere settling freely in a quiet medium to its radius 
s 

(r) and is given by the expression 


v 

s 


(p - p ) g r 2 

9 77 p m & 


in 


( 1 ) 


where and p^ represent the specific gravity of the sphere and the medium, 

respectively, g is the acceleration due to gravity, and ij m is the viscosity of 

the medium, v^ is also referred to as the Stokes’ velocity. The above 

expression for alumina particles (p = 3.965) settling in air (p = 1.22 x lo -3 

p m 

and ?? = 1.818 x 10 -4 poise) reduces to 

m 


v^ = 4.837 x 10 -2 r 2 (cm/sec) 


(2a) 


1. Adarsh Deepak, R. Ozarski, andJ. A, Thomson: A Standard Particle 
Sizer-Velocimeter. NASA Technical Memorandum, March 22, 1974 (to be 
published) . 

2. This effort was supported under NASA -Marshall Space Flight Center 
contract No. NAS8-24810. 



where r is in micrometers. For settling in water (p = 1.0 and7) m = 0.01 
poise), equation (2a) reduces to 


v = 6.463 x io " 4 r 2 (cm/sec) . (2b) 

s 


The numerical values of v versus r are given in Table 1 and are represented 
graphically in Figure 1. 

Aerosol particles accelerate to this limiting (terminal) velocity rapidly. 
The distance and time a spherical particle falls to reach 67 percent of the 
terminal velocity are as follows: for alumina particles settling in air. 


v 

Time = — = 0.493 x 10 “ 4 r 2 (sec) , 
g 


and 


Distance = — — 

2g 


1.192 x 10 -6 r 4 (cm) 


(3a) 


(3b) 


for alumina particles settling in water. 


v 

Time = — = 0.658 x lo -6 r 2 (sec) 
g 


(4a) 


and 


Distance = 

2g 


2.129 x io~ 10 r 4 (cm) 


(4b) 


Thus it takes only a few milliseconds and a fraction of a centimeter for the 
alumina particles (p = 3.965) to reach their terminal velocity. 
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TABLE 1. VALUES OF STOKES’ VELOCITY FOR ALUMINA 
PARTICLES OF DIFFERENT RADII SETTLING IN 
AIR AND WATER 


For alumina particles, 


Pp = 3.965 (in water at 25°C). 


For air, 


P m = 1.22 x 10' 3 


and 


I7 m = 1.818 x 1CT 4 poise (dyne-sec/cm 2 ). 
For water, 

“m = 1 ° 
and 

T? m = 0,01 poise. 

Stokes’ velocity in air, 

v g = 4.554 x 10" 1 r 2 (mm/ sec) . 

Stokes’ velocity in water, 

v = 6.457 x 10' 3 r 2 (mm/ sec). 


Particles 

of 

Substance 

p P 

Radius, 

r 

(pm) 

Diameter, 

2r 

(pm) 

Stokes’ velocity, 
v s 

Air 

(mm/sec) 

Water 

(mm/sec) 

AI2 O3 

3.965 

0.15 

0.3 

1.02 x 10‘ 2 

1.45 x lO" 4 



0.5 

1.0 

1.14 x 10' 1 

1.61 x 10' 3 



1.0 

2.0 

4.55 x 10’ 1 

6.46 x 10' 3 



1.5 

3.0 

1.02 

1.452 x 10' 2 



2.0 

4.0 

1.82 

2.583 x 10' 2 



2.5 

5.0 

2.84 

4.03 x 10' 2 



3.0 

6.0 

4.098 

5.811 x 10‘ 2 



4.0 

8.0 

7.286 

1.033 x 10' 1 



5.0 

10.0 

11.38 

1.61 x 10' 1 
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GENERAL DESCRIPTION OF THE METHOD 


Aerosol particles are allowed to fall freely in a vertical glass-walled 
cell (Fig. 2) in which the convection currents have been reduced to a minimum. 
The settling velocity of the particles is determined by photographing them while 
the light entering the lens is chopped at a known rate. The particles are 
illuminated by a 1- or 2-ram thick vertical sheet or slab of light projected into 
the cell. A camera aimed perpendicular to the slab of light photographs the 
particle tracks; typical exposure times are 0. 5 or 1 sec. The image of the 
falling particle is thus a series of dashes (with the chopper) (Fig. 3). The 
velocity of fall can then be measured from the spacing of the dashes. 

Under ordinary conditions the velocity thus measured is the sum of the 
Stokes’ velocity and the convection velocity. It is the experimental suppression 
of the convection currents and the theoretical elimination of their effects that 
is the major cause of the complexity of the experiment and of the limit of the 
applicability of the method to large particles. The procedures for damping the 
convection currents are described later. One finds it very difficult to reduce 
the convection velocity below about 0.5 mm/sec. However, it is possible to 
eliminate by subtraction most of the effects of convection by ensuring the 
presence of many submicron particles ( such as DOP particles approximately 
0.3 pm in diameter) . These move essentially with the convection air currents 
and can be used to measure the vertical component of convection velocities, 
which can then be subtracted from the vertical component of the total velocity 
of large particles to obtain their true Stokes’ velocity. 


BRIEF DESCRIPTION OF THE INSTRUMENTATION 


The instrument setup essentially consisted of a vertical glass-walled 
cell and a photographic system (Fig. 4) . 

The cell was a vertical glass-walled chamber which had airtight joints 
and can be closed-off airtight at the top and bottom. Its cross section can be 
octagonal, rectangular, or square. The one used for this research was 
octagonal (2-in. sides) in cross section and 18 in. high. Its walls were covered 
on the inside with (black) flock paper except for the appropriate 2-in. diameter 
openings for the passage of light and for the camera lens. This helped minimize 
the extraneous light entering the lens and provided a dark backdrop for the 
illuminated particles. The bottom was closed off by inserting a 1-in. thick 
aluminum plate. 
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(c) TOP VIEW 

Figure 2. Schematic illustration of the method. 
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Figure 3a. A typical photograph of settling AlgOg particles. 



35 mm SLIDE 


Figure 3b. Typical tracks for 0.5-sec exposure (schematic) . 
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Figure 4a. Photograph of the experimental setup. 

The plate can be cooled to and maintained at a steady temperature, the 
temperature of dry ice (-12°C) , over its entire surface by circulating a C0 2 
and acetone mixture through a passage inside the plate using a dry ice-and- 
acetone pump (Fig. 5). 

The photographic system consisted of a light source, a camera with a 
bellows, and a uniformly rotating chopper. The entire system was mounted 
between two parallel aluminum sheets (0.125 in. thick) held 5 in. apart by 
aluminum rods to form a sturdy unit for field use. 

The source of light was a 500-W slide projector (Bell and Howell) with 
a rectangular aperture (l mm by 1.5 in.) inserted in the position of a slide. 

The beam of light was focused by a short-focused convex lens (f/2) to form a 
sharp focus of the aperture at the center of the cell. Near the focus a vertical 
sheet or slab of light 1 or 2 mm thick illuminated the falling particles. The 
depth of focus was about 5 mm. 

The camera unit was an f/1.7 Minolta lens mounted in front of a bellows, 
and the unit was fixed to the external framework (Fig. 4) . The bellows was 
adjusted to obtain a sharp image of an illuminated object at the position of the 
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F igure 4b. Schematic diagram of the instrumentation. 



projector beam focus, approximately 3.5 in, in front of the lens, so that maxi- 
mum magnification could be obtained. The magnification in this case was 
approximately 0.8. The camera body could be removed from the bellows unit 
for film changing. 
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A uniformly rotating chopper was used in front of the lens to accurately 
select and measure the particle track lengths in the photograph. With the aid 
of a belt-and-pulley arrangement the chopper (2 or 4 blades) was rotated by a 
small synchronous electric motor (300 rpm) fixed to the external framework. 
By choosing an appropriate combination of pulleys and chopper blades, one 
could obtain different chopping rates (typical values were 10 or 100 chops/sec). 

Tri-X (ASA 400) film was- found to be the most suitable film considering 
speed, resolving power, and the ease of developing. It is sensitive (fast) 
enough to record tracks of particles of diameters as low as 0.3 )im often moving 
with velocities as high as 3 mm/sec, when the bellows was set for a magnifica- 
tion of approximately 0.8. 


A SETTLING EXPERIMENT WITH ALUMINA 
PARTICLES 


An extensive series of experiments using alumina particles was carried 
out to test the performance of the method and to define the limits of its accuracy 
and applicability. The samples of AlgOg particles used were not strictly 
monodisperse samples, but had size distributions which were narrowly poly- 
dispersed about the mode diameter values of 0.3 Mm, 1.0 Mm, and 3. 0 Mm . 3 
These samples were chosen over the more accurately calibrated particles 
(such as latex particles) because they were relatively inexpensive and provided 
enough calibration for the purposes of making the initial feasibility and per- 
formance studies of the method. Only one set of results of a typical experi- 
ment with 3.0-Mni particles is presented in this report. 

Aerosols of AlgOg were generated by impinging a jet of air on a small 
sample of AlgOg particles placed inside a long vertical tube (3-in. diameter, 

5-ft length) . The AlgOg aerosols were drawn off into the glass cell as needed. 

After the particles were introduced into the cell, the top was closed 
airtight. After turning the projector light on, several bright specks of light 
could be seen through the camera viewfinder swirling around as they traversed 
through the beam. These random flow patterns were created to some extent by 
the initial disturbance created while introducing the particles in the cell, and 
to a large extent by the convection air currents. These convection currents 


3. The data were provided by the supplier, Adolf Meller Co. , Rhode Island. 


10 



arise because of the small differences in the temperatures of the various parts 
of the cell. As mentioned before, the velocity of the particle is then the sum 
of the Stokes 1 velocity and the convection velocity. 

It is essential to the effective working of the method to suppress con- 
vection air currents. One resorts to several procedures to reduce the currents. 
The simplest procedure is to make all the joints of the cell truly airtight. The 
best test for air leakage is to check for water leakage. Also, closing the top 
of the cell helps. 

The most effective procedure for damping the convection currents is to 
cool the bottom of the cell with ice, dry ice, or dry ice and acetone. This 
creates a stabilizing temperature gradient, and consequently density stratifica- 
tion, in the column of air in the cell. The best results were obtained by using 
the acetone-dry ice pump to cool the plate to -12°C. The motion of the falling 
particles was seen to quickly stabilize to a nearly vertical direction of fall 
within 90 sec after the insertion of particles. 


Inspite of taking all these precautions, slow motion convection currents 
(~ 1 mm/sec) still seemed to persist. One can measure their velocities by 
introducing 0,3-pm diameter DOP particles (Royeo DOP generator) into the 
cell. These particles essentially follow the convection currents, since their 
Stokes 1 velocities (<*' 0.005 mm/sec) are negligibly small compared to the 
convection velocities (typically 1 mm/sec). The same is true with 1.0-pm 
or 0.3-pm diameter AljOj particles. It is difficult to distinguish their tracks 
from those of DOP particles, making the measurements unreliable. The DOP 
particles have been seen to follow the currents in the airtight cell for periods 
of over 8 hr. The tracks of 3. 0-pm diameter A^Og particles, however, are 
easily distinguishable from those of DOP aerosols. It appears that had we 
been able to use 2.0-pm diameter A1 2 0 3 particles, then, with some care, their 
tracks could also be distinguishable from DOP tracks. Thus it appears that 
the lower limit of diameters for which this method is accurate is 3.0 pm, but 
the limit can be, perhaps, decreased to 2. 0 pm for particles whose specific 


gravity is approximately 4.0 (or 



for specific gravity p ) . 


When the particles stabilized to a nearly vertical direction of fall, 
photographs were taken for 0.5-sec or 1-sec exposure times. Several photo- 
graphs for each sample were necessary to ensure statistical accuracy of the 
results. 


It is important to periodically check during the experiment and make 
certain portions of the glass and the camera lens do not become fogged. 
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THE PHOTOGRAPHIC DATA REDUCTION 


The negatives were projected upon a screen and the lengths of the 
selected particle tracks were measured by hand. The visual distinction between 
the DOP tracks and the A^Og tracks was made by noting that all of the DOP 
tracks in a given region of the frame were of more or less uniform lengths and 
had the same inclination to the vertical, while the tracks of the heavier AljOg 
particles were larger and less inclined to the vertical. Notice that the hori- 
zontal components for all the tracks in the given region should be the same 
length. An exception to this would be a very large particle (typical diameters 
> 20 . 0 pm) which might travel almost vertically through the field of view. 

The selection of the AlgOg tracks suitable for measurements was made 
as follows. (A few tracks that appear typical are schematically drawn in 
Figure 3b.) For a chopper speed of 10 chops/sec, one expects, for a 0.5-sec 
exposure, a series of five dashes to represent the track of each particle if the 
particle falls vertically and remains within the 1 -mm thick illumination slab 
during the 0.5 sec. If however, the particle trajectory is transversely inclined 
to this slab of light such that the particle remains illuminated for only a fraction 
of the exposure time, then it may leave a track such as P 2 (Fig. 3b) having 
a number of dashes less than five. If the number of dashes is less than three 
(for 0.5 -sec exposure) the particles are excluded from measurement. How- 
ever, it was found that under the best working conditions, the fraction of such 
particles is rather small and their exclusion introduces a small error in the 
particle number density results. Thus the chopping of light not only helps one 
make an accurate measurement of particle velocities but also aids in the 
selection of particles that remain illuminated throughout the exposure. 

Ideally, the convection velocity component should be measured in the 
vicinity of each of the different Al^Oj tracks by selecting two DOP tracks 
preferably close to and on opposite sides of an AI 2 O 3 track. This is illustrated 
in Figure 3 b by the DOP tracks, D t and D 2 , about the Al^Og track, P 0 . The 
lengths of the horizontal components (d t , d 2 , d 0 ) of the three tracks should be 
equal. Then the true Stokes’ velocity component is 

Oh*-) 
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where h 0 , hj, and hg are vertical components of their track lengths. The 
Stokes' velocity 


v = h /t 
s s 


where t is the exposure time. Using the Stokes’ law, one obtains the radius, 

r. 


DATA ANALYSIS 


To analyze the data, the particle tracks in a number of photographs are 
measured and reduced to particle velocity (and consequently particle radius). 
To obtain a 10-percent error bar in each of the 10, for example, particle 
radius intervals, it is necessary to reduce about 1000 tracks (± 20 percent in 
10 intervals requires about 250 tracks) . 

If the particle tracks are all entirely contained within the viewed 
volume, the average particle size distribution can be evaluated as shown below. 

The height of the viewed region = h. 

The area of the horizontal cross section viewed region = A . 


The volume of the viewed region = hA . 


If Np (r^, Ar) is the number of particles seen in the pth photograph 

A r A r 

for which the particle radius lies between r. - -r- and r. + —— and ]SL is 

i 2 l 2 f 

the total number of frames viewed, then the size distribution for the ith radius 
interval is given by 


= rl l 

f P 


N (r., Ar) 
P i 

hA Ar 
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fn(r)]. N ua Ar ^ N p ^i’ ’ 

f p 

if h. A, and Ar are constant. 

For particles whose tracks do not both originate and/or terminate in 
the illuminated volume, the procedure is to count the number of tracks, 

N (r , Ar; z ) , which pass through a given elevation z. in the cell during 
P i j 3 

the exposure time. If measurements are made at the number (N ) of levels 

in the photograph, then the particle size distribution in this case is given by 4 



where v (r.) is the velocity and At is the exposure time. N (r., Ar;z.) 

s i P P i J 

is the number of particles crossing lines at each of the z^ levels, regardless 
of whether the same particle has been counted on other levels. 


DISCUSSION OF RESULTS 

The results of experiments with only the 3, 0-pm diameter Al^C^ 
particles are presented here. The measurements of track lengths were done 
by hand from the negative which was enlarged by projecting on a screen. The 
number density of the particles in each of the photographs analyzed is about 
175 cm -3 . The size distribution is given in Table 2 and Figure 6. The peak 
of the n(r) versus r curve is at about 1,45 pm. This agrees well with the 
value of 1.5 pm supplied by the manufacturer for the sample of alumina 
particles, which are supposed to be narrowly polydispersed about 1.5 pm. 

The results of experiments with A1 2 0 3 particles of diameters 0.3 pm 
and 1.0 pm were inconclusive because it was found difficult to distinguish the 
AlgC^ tracks from the DOP tracks in the photograph. However it is felt that 

4 0 Personal communication with Dr. J. A. Thomson, Physical Dynamics, 
Inc. , Berkeley, California. 
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TABLE 2. SIZE DISTRIBUTION OF SETTLING ALUMINA PARTICLES 


Volume of the illuminated particles is 0.04 cm” 3 . 

Average number density of the AljOs particles is 175 cm -3 . 

Average number density of the DOP particles is 5000 cm -3 . 

Radius Range 

(m) 

n(r) , cm -3 pm -1 

0 - 0.5 

0 

0.5 - 0.75 

14.3 

0.75 - 1.0 

28.6 

1.0 -1.25 

155.0 

1.25 - 1.5 

186 . 0 

1.5 -1.75 

72.1 

1.75 - 2.0 

85.0 

2.0 -2.5 

21.8 

2.5 -3.25 

9.28 

3.25 - 3.75 

9.46 

3.75 _ 4.25 

4.00 


for 1.0 pm particles, it would be possible to get meaningful reliability by using 
more sophisticated optical track measurement equipment. 

Even though one can photograph 0.3-pm diameter DOP particles moving 
with speeds of approximately 1.0 mm/sec, it is thus the speed of the convection 
currents that defines the lower limit of the diameters of aerosols for which this 
method is accurate. For particles of alumina, it seems possible to push the 
lower limit of diameters to 2. 0 pm, and for particles of specific gravity p 

P 

to 2 (v4 7p ) pm, without loss of accuracy. 
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